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(57) The present invention is directed to a nucleic 
acid quantification kit and method for determining the 
initial concentration or mass fraction of a target nucleic 
acid present in a sample. Illustrative embodiments in- 
clude real-time competitive quantitative polymerase 
chain reaction (PCR) to determine the copy number or 



mass fraction of a target nucleic acid sequence in a sam- 
ple and use of a thermodynamically based signal 
processing algorithm, with or without PCR, to provide 
mass fraction information. 
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Description 
BACKGROUND 



5 fOOOIl The polymerase chain reaction (PCR) is a technique of synthesizing large quant.t.es of a prese ected °NA 
segment The technique 

faboTatorv PCms achieved by separating the DN A into its two complementary strands, bind.ng a pnmer to each single 
TaZ Xe end ol the given DNA segment where synthesis will start, and adding a DNA polymerase to synthe «e 
Secomplementarystrando 

10 number of copies of the selected DNA segment have been synthesized. , h ^ tom 
Zofl During a typical PCR reaction, double stranded DNA is separated into its single strands by raising the tem- 
perafure of the 5nA containing sample to a denaturing temperature where the two DNA strands separate (i.e., the 
C n temperature "the DNA") and then the sample is cooled to a lower temperature that allows the specfc primers 
to aiach SaT) and replication to occur (extend). In illustrated embodiments, a thermostable polymerase ,s utilized 
S^SZ252n reaction, such as Taq DNA Polymerase and derivatives thereof, including , the , Stoffe. fragment 
of Taq DNA polymerase and KlenTaql polymerase (a 5 --exonuclease deficient variant of Taq polymerase - see U.S. 

SolTTne^ 

One of the mafor reasons for this increased use of quantitative PCR derives from the fact that PCR has • 
?ve ordei of magnitude better than the best blotting procedures. This sensitivity makes PCR as a ; q^ntttative too. 
h ghly desirable, however, the use of a system undergoing exponential amplification ,s not <^£T 
tion. Small differences between sample sizes can become huge difference in results when they are amp.if .ed through 
forty doublings. , 
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25 Kinetic PCR 



T00041 Atypical PCR reaction profile can be thought of has having three segments: an early lag phase, an exponential 
growtlph^ 

signal S the probe system used to detect the PCR product. The exponential growth phase begins when sufficient 
30 product has accumulated to be detected by the instrument. During this 'log" phase the ampl ication course , . descr bed 
bv the eouation T = T (E) n where T n is the amount of target sequence at cycle n, T 0 is the initial amount of target, 
and E is the Sidency of amplification. Finally, in the plateau phase, the amplification efficiency drops o« extreme y 
rapidly Product competes more and more effectively with primers for annealing and the amount of enzyme becomes 
limitino The exoonential equation no longer holds in the plateau phase. 
3* 00 ] Most o'f the quantitative information is found in the exponential cycles, but the exponential eye es > 

comprise only 4 or 5 cycles out of 40. With traditional PCR methods, finding these informative cycles requires ^ that the 
Son be split into multiple reaction tubes that are assayed for PCR product after varying numbers .of 
requires either assaying many tubes, or a fairly good idea of the answer before the experiment ,s begun Once the 
posin of the exponential phase is determined, the experimental phase.can be compared to known standards and 
40 the copy number can be calculated. 

Competitive Quantitative PCR 

[0006] Competitive quantitative PCR methods were developed to attempt to overcome diffteulties ^ s te SO r C Jn2m«ed 
indinq the exponential phase of the reaction and to obtain greater precision. A competitor sequence rs constructed 
hm s am^imed using the same primers as are used to amplify the target sequence. Competitor and target are differ- 
em a ed usually by length or interna, sequence, and the relative amount of competitor and target are measure i after 
amplLtion . If L target and the competitor are amplified with equal efficiency, then their ratio at the end o the reaction 
w7be the same as the ratio had been at the beginning. This holds true even into the plateau P»^£ng«£Ji 
decline in efficiency at the same rate. Thus, finding the exponential region is no longer a prob ^iT^r e^ 
in the same tubes with the unknown targets allows for additional control not possible with kinetic methods^ For example. 
Lading the competitor before mRNA purification would control for variations in sample preparation and reverse tran- 

[o'ooT The use of currently available competitive PCR techniques continues to suffer from several di^FHI- 
y . the competitor sequence must be constructed to be as similar as possible to the target seque ,nc, ■ ° h * 

efficiency of amplification, yet the two sequences must be distinguishable from one another. If the compel or is too 
ST^^c. to the target, heteroduplexes form during the PCR that skew the ratio of the product to ta^k. 
[0008 Tn addition, competitor must be added to the unknown sample at a concentration approximating that of the 
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■■Hi 

problem of forming heteroduplexcs during the PCR react.on. 
Real-Time Quantitative PCR 

[0010] Developments in instrumentation have now made real-time monitoring of PCR reactions possible and thus 
\ZT™ - 0. monitoring a. o.oh oyo.e ,or qua «a«v. PCR was ^>°*><*> H|« i 0.^ 

EST '^e'sa SslZanSn'anallu^scanl m on, t o* g ,»hni qU s S r,av,^ak,na,ic PC R »,gnir« ! ,n 1 ,,,a,ier,han 



invention. 

SUMMARY OF THE INVENTION 
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and/or melting curve analysis. nnmnf >x\tor nucleic acid sequence is prepared to have the identical sequence 

[0017] in an illustrated *^ dme ^*^ , oca ted at an internal position on the com- 

as the target nucleic acid sequence w, h the ««pton ^^^™^ 6b compos ition as the corresponding 

^-"Othe.en—.^^ 

in length an*., composition, but the competitor and target 3^Xt31«ne«pe«nent/tlon. 
same efflcienc,. Such amplification etliciencea can be °" = G ,^'" me competitor nocleio acid 

nucleic acid sequence, and the target probe .s con ''9 ured »° ^'"^ Accordingly when the anchor, target and 
corresponding to the unique region of the competitor ^^^^^^^ and competitor nuc.eic 
competitor probes hybridize to their respect.ve complementary ZlhS Ts well as the donor fluorophore and the 
acid sequences, the donor f.uorophore and the f.rst ^Ptor f.u^opho " the me P asurem ent of 

second acceptor f.uorophore are placed ,n a resonance «^y tra" star . ^^^^m^ of the target nucleic 

rc^r.r,^^ 

i^^tz^zz^^r^v^ r,:^ ~ .he 
^j^zx^zs^szzz* - - - — — * **** 

for target and competitor, logC 0 = logE(An) + log i c , wnere o amount of target. Because this 

rheT.~^le a ^ 

melting temperature, illuminating the test sample monitor ng fraction of the target nucleic 

using a thermodynamic^, based signal P'^'^.^^^^^^JST^ „ different from the 
acids. The internal standard may consist of an artifiaa. compear 

target nucleic acid sequence by oheor more basesjf '^^^^'^^^„^ ag 0 i M 
then the initial copy number of the target nucleic acid s ^"^^^ l ^^^ ummMion ma y be in deter- 

x^:e°ir:rrpt^ 

sected tissue samples from cancer patients. ™ QrW t n those skilled in the art upon consideration 

SLoC-lS^ 
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as presently perceived. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0025] 

M two «-W P«*«« «> «*><*»' »' a 

acid seauence resulting in fluorescent resonance energy transfer; „ Mrtro<s 

S 3a b represent typical external standard curves using hybridization data. F.g. 3a is a plot of the log fluores- 
Figs. 3a b represent wroi e vs second derivat)Ve maximvjm; 

FrVatre^eL 

S^T^rl pTotted against temperature (Fig. 4b), with homozygous mutant (.•••), homozygous w,.d type 

fcsla-^ Uic acids. Fig. 5a shows melting peaKs generated from a 

Sint our^J The area under each curve is calculated using non-linear regression to fit the melt.ng peak data to 
a Ss" TcurTe 5b shows v ari ous amplification curves on a log fluorescence vs. cycle number p.ot. Fjjfc 
shows the data from Fig. 5b converted into a ,ogC 0 = log E (An) + logT 0 curve (solid line show cross.ng pomts from 

^ "e^r^^^ sequences of the HER-2/neu (target), its competitor, the reporter and anchor 

pmbL The "rented Siting tenures Tm of the reporter probe hybridized «o either the target or compet.tor 

Rn ^"diagrammatical representation of the strategy used to create the competitive DNA fragment for HPV 16; 
Fig 9 ts a diagrTmmatica. representation of the hybridization probes used to detect the internal quant,f,cat.on 

STSE a qraphTc Z^on^SS^n efficiency of the Interna. Quantification Standard (A)and Artificial 
HPV ' 16 temple ^ (pThe dlta are presented as the average «f at least three separate data points, w.th standard 

R^Ta^lSe atypical interna, control reaction demonstrating fluorescent dat^ 
P^aTareaS^ 

of 1x104 The internal quantification standard is at mrtial template concentrat.ons of 1x10 (F.g. iza), I 8 

Sou; Xentnfre^iondata points. The 95% confidence interva, at each ratio of competitor to target ,s 

fT ™ between me.ting peak area and product concentration for mutant and wild-type 

HER iS^SL^SLi by hybridization probes using me.ting curve analysis software. Art,f,c.a. otago nudeot.de 
template: werTrSxed wih probes at various concentrations and me.ting peak-area was deferred us.ng L.ght- 

curves were analyzed by the TMBSP algorithm to determine the rat.os of mutant and w.ld-type PGR products 
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the FRET pair probes (Figs. 1 6c. and d); Fig 1 6a and c show melting data (fluorescence vs temperature) and Figs. 
17l>andd8lK.v.lhemeningpeakdala(r.egalive«tslderlvative.clF/dT). 

approximation algorithm is based on this model to approximate the melt.ng curve. 
15 DETAILED DESCRIPTION OF THE INVENTION 

to be quantitated using standard techniques. The method uses a com ^ e J ns o( calcula ting the initial 

amplification or melting curve ^£Z££ZZZZ^ ^n limited* quantifi- 

concentration of the target sequence. To date, all real-time POH quant, ca hw ube -to-tube differences 

cation relative to an externa, standa^^ 

the sample container 450, .s controlled w.th ^^^ ^ ™ ( J drjven via a shaft 496 and a motor 494. 
fan 498, which is provided to move the air ,n the Seated P ath 472 J h ^ n ' frustrated embodiment, twenty-four 
The fan 498 forces air into the aperture 470 and f« the heating cartridge 

sample containers 450 (two of which are represented .n F.gJ) ^ symmet ncahy ar rang ^ ^ 

474 and the central duct 476. The sample containers 450 are ™ " connected to the motor 488 

carousel 480 is positioned by a stepper motor 488 provided w. * J^^^JJJ w h ichtc.udes ah excitation 
via ashaft486. Fluorescence from each sample container ,s 0 ' '^^^ .™ , ^ ound in U.S. Patent 
radiation source 468 and photodetectors 464 and 466. M ° re othertherma. 
Application No 08/869,275. It is understood that this described embodiment is merely exemplary ana 

at least as low as the first temperature for a pre-determ.ned length of T^^^T^Z^ number of times, 
then raised back to the second temperature, and the bio.og.ca sample ^^^^^JJSii. profile wherein 
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1997). These detection (oimats include but are not limited to: 

1 Use of fluorescent double-stranded DNA recognizing compounds (see Figs. 2a-b) 

2. Taq Man principle (see Figs. 2c-d) 

of accumulating amplification product. 
3. Molecular beacons 

See, for example. U.S. Patent No. 5,118,801. 

4 Increased FRET upon hybridization (see Figs. 2e-f) 

mmmwmMm. 

'ZSS^S^iSSSSZ — n 0, ft. dono, ftoiety and ™as u renne n , o, nuores- 

mentary to the sequence o the anatyle. to the ana,yte under appropriate hybridization 
of resonance energy transfer is 



(8.785E 5 ) (f 1 ) (k 2 ) (n* 4 ) (q D ) (R 6 ) (J D A>. 
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t= excited state lifetime of the donor in the absence of the acceptor; 

k*= an orientation factor between the donor and acceptor; 

n= refractive index of the visible light in the intervening medium; 

q n = quantum efficiency of the donor in the absence of the acceptor; 

R- distance between the donor and acceptor measured in Angstroms; and _ 

^JS^^Fo) <e A ) M with respect to W at all overlapping wavelengths w.th: 

F D = peak normalized fluorescence spectrum of the donor; 

e A = molar absorption coefficient of the acceptor (M^cnr 1 ); and 

W 4 = wavelength (nm). 

,004,, For an, given don., nnd accept., a dis,ancn »ha,a ^^T^^e^T,!^^ 
and * ab.,avia,ad R 0 . Bacaose ,h. ,a,a o, resonance ane irg, »en*er dec** £ *> .a* pow. 

=Ta£r^^ 

between the donor and acceptor fluorophores is about 1 nucleot.de. rpr obe-primer"), then the two 

I0 043] When one of the .abe.ed oligonucleotides also * n «™*"*™%£" ^donor and acceptor, fluoro- 

a detectable increased fluorescence by the acceptor fluorophore. 

~zzz=z=zz^: ™=- 

5. Single-Labeled Oligonucleotides 

PM, S,n g ,e-,abe,ed o„onuc r „das a,a ^'^1^2^ 
MlgonnoWlldaa may be osed indapendenlly ol any """"f' Appfcaoon No. 09«27.842. mad 

Sr. «,and, hvbndb.don may reap,, in aim., a decease o, increase ,„ lluo.eaeence. 
Probe Systems for the LightCycler® 
(0048, Aseguencespecificprobesysternfor^^^ 

wherein two fluorophores of a FRET pair are brought close ^^^^^^ between the PGR 
energy transfer occurs (see Figs. 2e-f). Two adjacent hybnd.zat.on T^ 8 ™^"*™,* Z .acceptor fluorophores. 
pnme.i,one.abe l edatthe3'endwithadonorf W ^ 
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donor dve fluorescein and two different acceptor dyes, LightCycler Red 640 (LCRed 640) and LightCycler Red 705 
Sea %TwMeF*ET oligonucleotide pairs are commonly used with the LightCycler® and are used vanous ex- 
Lpls herelnTis understood that other sequence specific probes may be used within the scope of th.s .nvenfon. 

Real-Time Kinetic PCR on the LightCycler® 

rO0491 The LightCycler® can be used with either double stranded DNA binding dyes such as SYBR™ Green I or 
ybr dizalon presto monitor the PGR reaction. Fig. 3a and 3b show typical ^.standard curves . us,ng h*md- 
iTation orobes The donor probe was labeled with fluorescein and the acceptor w.th LCRed 640. The data are plotted 
as Ihe ra°o of acceptor to donor fluorescence. The initial concentration of standard ranged from 1 0* to 1 d cop.es of 
target per 10 u.l reaction. 

Mutation Detection using the LightCyclei® 

ro050] Monitoring once each cycle provides useful information for quantification. Additional information ^ available 
f ZLcence is monitored continuously during temperature transitions. The hybridization state of the probes can be 
iSZZ^me^STS- change inLrescence as the temperature is varied. Hybridization probe melt.ng occurs 
at a characteristic temperature that can be exploited for product identification and mutation detect.on. 

Quantification by Kinetic PCR 

[00511 The temperature dependence of the fluorescence from hybridization of the probes may be demonstrated with 
Celcence vs temperature plots (Fig. 4a). The illustrated plots were generated by monitoring a s.ng.e sample every 
oT^during a slow (0.2-C/second temperature ramp from 45°C to 75°C. The product is denatured and hen rapidly 
cooled ( ucc/second) to 45o C . At low temperature the probes hybrid^ 

ratio (for example LCRed 640/ fluorescein) increases. During heating, the probes d.ssoc.ate n the 55 to 65 C range, 
ret- n florescence ratio to background levels. The derivative of this curve is -'-' a ;^;; t h ; e ^ e o C 
ature and plotted against temperature (Fig. 4b). This produces a melt.ng peak centered around the T m of the probe. 
Discrimination based on hybridization temperatures is a powerful tool for mutation detect.on. 

A Method Combining Mutation Detection with Quantification 

[0052] The use of an interna, standard in competitive quantitative PCR.assays involves careful selection oft the com- 
'oetito used as the internal standard. The competitor and the target in compete quantitative PCR assays must MhU 
con radictor^ criteria The two nucleic acid sequences must amplify with the same efficiency, generally requmng them 
to be a^ Slar as possible. But they must alo be differentiab.e and not prone to heterodup.ex format.on. requmng 

ISu^ ^he 8 utrnate in similarity between target and template is a single base pair change. 
hat a sinole base change would have a significant effect on efficiency of ampl.f.cat.on. In accordance with one em 
Z iment o "tEK£S£ the LightCyc.e^is used to differentiate between a target and a competitor diff ering , by onj 
a sS base pair, as in a single base'pair mutation. Under proper condition's, hybridization probes detect only one of 
the DNA strands, so heteroduplex formation during amplification does not affect the results. 

0054] fn the course of the development of the LightCyclei®, software has been developed for analysis of reaM me 
luorescence data Fig. 5a is a representative melting curve. The software calculates the area under each curve using 
non-Hnear regrei'lon^o fit the melting peak data to a Gaussian curve. This module serves as the basis ; of tN , so«wa e 
to^ZJL* using the T m method. The relative peak areas of target and compeftor are used to calculate the 
relative amounts of the two products. nn . r . For each curtfe tne 

[0055] Fig. 5b shows various amplification curves on a log fluorescence vs. cycle number plot For each cu-ve _the 
point in the amplification curve where the second derivative is at a maximum is identified, that is the po nt of maxima, 
fn rease in thTra te of increase. This f ractiona. cyc.enumber is used to describe the position of the amplification cwve. 
Un fke tradi ionarthreshold" methods that define the curve position relative to background no.se this approach allows 
mc au toma I determination of the positions of the amplification curves based on the shape of the curve. See U.S. 
Patent No 6 38 module serves as the basis of the software for the mu.»i^o.ormethod The re.; = J 
oftrget and competitor are determined by .coking at the fractiona. cyc.e difference in the pos.t.ons of the two amp..- 
fication curves, as shown in Fig. 5c. 
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A Method Combining Kinetic PCR with Internal Standards 

rnnsfii m an alternative embodiment, the competitor/internal standard is distinguished from the target nucleic acid 

a "super bright" blue light emitting diode. Fluorescence signals are acqu.red from photod.odes after bandpass f.ltermg 

, 5 iooBr™ - * concern, as on,y one o, the DMA strands is detected by 

hiridizat on proves Work with external standards has shown that the position of ampl.hcat.on curves .s more 
the hybr.d.za ion P^ S g W °7J 1 , eve , s Accordin g ly , since data are collected every cycle in this kmet.c method, 
memor TrSZ da a ^ ~Sled. Advantageous*, the present method does no, depend on a single 
r e a"ule« detl the product Los. Instead, the relative positions of the entire amplification curves are used to 

20 ZsTl SlllS^e target and competitor concentrations give amplification curves that are in the 
l ie , rnSSm in the curve position between target and competitor can be usee I to calculate the rat.o of 
fargc. and competitor. This method provides precise estimates of target and competitor amounts. 

25 Delia C.T. Equation Determination 

roOSOl The above approach has not previously been used with quantification with interna, standards. Thus a con- 
lentel, preLably linear mathematical re.ationship between the target and the competes curve posrt.ons _and the* 
. ".Xe concentrations is needed. If target and competitor have the same efficiency, then at the second denvat.ve 
30 maximum for the target: 



T n « = T 0 (E) nt 



35 



40 



45 



where T . is the amount of target at the second derivative maximum, T 0 is the initial amount of target, E is the average 
eSnc?, of theTeacL, and nt is the fractions, cycle number of the second derivative max.mum. S.mHariy at the 



second derivative maximum for the competitor: 

C nc = C o( E) 



nc 



where C is the amount of competitor at the second derivative maximum. C 0 is the initial amount of competrtor, E is 
^ «em« ^1^2 the reaction, and nc is the fractional cyc.e number of the second derivahve maximum. 
00611 The ^sSond derivative method is sensitive to the shape of the amplification curve, not the absolute fluores- 
cing }^ J ^£^TL amplification curve should not be significantly affected by differences^ „ , .gna.ng 
cence level, me pa I Red 740 The point where the second derivative is at a maximum does not reflect 

maxima, the concentrations of target and competitor should be equal. Therefore: 



50 C nt =T nc 



And so it follows that: 
55 ' . C 0 (E) nc =T 0 (E) nt 



Rearranging: 
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cyr=(E) nt / ( E) nc 



log(C 0 H- 0 ) = log[(E) nt /(E) nc ] 



logC Q - logT c = ntlogE - nclogE 
logC Q - logT D = logE(nt-nc) 
nt-nc is the cycle shift between target and competitor which we can call An, substituting: 

logC Q - logT c = logE(An) 

And rearranging: 

logC c = logE(An) + logT 0 
u f^r^ <i my 4- b so a olot of the initial competitor concentration versus the 

to the log of the initial target concentration. 



EXAMPLE 1 

(00631 The totlowing e*per,men, is conducted to ddntirm ,hat eaua, ooddentratiods ot initial targe, add competitor 

are the same tor equal concentrations. This •^h^STSSSS^I^ S 4 IMhe dMerenoe is no, zero, 

difference between the two channels. 
EXAMPLE 2 

l0 065] The following experiment is conducted to confirm that the dynamic range of the assay is at least an order of 
magnitude on either side of the target concentration. concentrated product will reach a 

fCTminlmum JL* range o, a, leas, a "^-^J™^^ p. detected is tested, 
^eo^per^^^ 

and the precision of the calculated target number ^ de ^ ned established witn 10 * copies of target, the max- 
[0068] Once the maximum target to competitor d'fference has been ^^"^^J^ri Targe t fr0 m 1 0 1 to 



11 



. <EP 128831 4A2_L> 



EP 1 288 314 A2 

of the calculated target number is determined. 
EXAMPLE 3 



10 
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[0069] The following experiment is conducted to determine the effect of target copy number on the accuracy and 

nSSTS^TSm PGR experiments are analyzed for precision and accuracy. For each starting copy number of 
araet from 1^ to 106% 95 % confidence interval is calculated. The inter-assay and intra-assay precis.on ^o«J- 
cSe bTmeasuring the coefficient o, variation (%CV) within and ~P%^ 

~ f tomDl f rnm im to 1Q6 At 10^ or 102 copies, it is expected that the %CVs will be around 100 /<>. At the nigner copy 

be around 25%. A 25% CV would allow easy discrimination of two-fold differences. 



40 



45 



Software 



is r0071l The curve positions are calculated using the second derivative maximum method. This method, which de- 
pends on curve shape and not absolute signal, is believed to be more resistant to differences 
h!t«,P P n the channels The cycle shift is plotted againsUhe initial competitor concentration and a line .s fit to the data. 

'^S^X^^ -™ a " ie — s < %cv < 50) - lhen the spftware supports this ca,culat,on as w 

20 EXAMPLE 4 

100721 A method for real-time competitive quantitative PCR in the LightCycler® using a competitor which differs f rom 
EJget by onH single base is described in the following experiment. The target and the compear are d.stinguished 
by the differential melting of fluorescently labeled hybridization probes. 

Experimental Design 

100731 The target for quantification in this example is the human HER-2/neu gene. The HER-2/neu gene is amplified 
fn 25? of breast tumors and the degree of amplification (usually 2-50 fold) correlates with survival time Fig. 7 shows 
a oSg oTpTobe Tor HER 2/neu. With this design, the competitor has a CA mismatch with the hybridization probe_ 
A ?A m^matoh in the center of a probe results in aT m shift of 5-10°C, sufficient to allow for separation of the matched 

'he P-ers that flank these probes (not shown) were designed using the Primer 
Designer™ software (Scientific and Educational Software). 

35 Construction of the Competitor 

r00741 Wild type HER-2/neu PCR product generated from human genomic DNA is used as the target. The competitor 
s aeneraS by amplication of HER 2/neu from genomic DNA with a mutagenic PCR primer containing a G to A change 
2X^£7tESr Products are gel purified, diluted, and then reamplified with the amplification Pnmers^e 
J , Z Irl Z, n„ ifiPd and used as taraet and competitor. The introduction of the mutation is confirmed by sequencing. 
ffSCSS aredetermined by Molecular Probes PicoGreen dsDNA quantification 
assay or by the limiting dilution method as discussed above. 

Probe Synthesis and Purification 

roorei The probes are shown In Fig. 7. The anchor probe is 3' fluorescein labeled. The acceptor probe is labeled on 
[he 5 ?.nc ISmJ itCyder R" 640 and is blocked on the 9 end by a phosphate. Probes are synthesized and purif.ed 
as discussed above. 

so Quantification with Internal Standards 

r00771 First a determination is made that the signals from target and competitor (that is, the melting peak areas) are 
Smon^o the r o u n nt of target present. This is first done with purified PCR P-ducts. Wild type 
HER-2/neu are mixed in equal concentrations from 10«> to ^ copies per tube. The meltmg peaks are obtained by 
55 TlZrow ZgZ temperature below the annealing temperature of the probes and then slowly heating , (0.2 -a 
Lecon^) to a temperature 15'C above the melting temperature of the probes. Fluorescence is acquired every 0, C 
durmg X rl^ZLo of the areas under the best-fit Gaussians is compared to the known .nitia. target/competitor 
ratio of 1 .0. Statistical tests produce a ratio of 1 .0 that falls within the 95% conf.dence intervals. 
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ES3KS3S3S3SSSESS 

SanyUa es but will not exceed the amount of primer used. Hybridization probe "T^^^l^SS! 
n i ™ and o 5 CM orimers so the highest concentration of product that can theoretically be produced by PCR would 
be SL een 0 Z 5 5l P e.iminary experiments indicated that accurate measurement of product : amour*, by 
m^^Peak areas needed probe concentrations in excess of the total amount of PCR product produced after ampl,- 
SSThspSd problems for the standard LightCyc.eKB) optics, ^^^^S^TTt.^^ 
™S exceed the detection level in the F1 channel. To over come this problem the F1 optics of a L «MCyd«« 
was *£SZZ biock -90% of the f.uorescent signa. transmitted to the F1 detector. '"S^SS^JX 
tions of probe could be used so that the probe concentrations are always .n excess of product. Reconstructed I me It ng 
exoeriments using artificial templates of known concentration were designed to measure peak areas with this modrfied 
rnsSmSxiss probe Fig. 1 4 shows that there is a linear correlation between melting peak areas and produc 
LotceTra ions ^between 0 1 and 0 4 CM using 1 .0 CM of each probe. These results indicate that end-point PCR product 
^ptTc£ESo„ of 0.5 pM or .ess) will consistently produce melting peak areas within this ..near range and 
yield quantitative information. Dynamic range of quantification by melting peak a^lys.s ' . (or a ten . 

rooeoi A linear relationship between melting peak area and amount of PCR product could be established for a ten 

itS^XraUnq ana.ysis software. To broaden the dynamic range of this technique, a novel method of melting 
eJ^^^^toTi based on a Thermodynamic Modeling based Signal Processing (TMBSP, see Example 
8,7it?^2 *. components of a heterogeneous melting curve are quantitatively described ,n terms 
J thAir volume fractions with respect to homogeneous melting curves for each component. 

?0u81 Fiqure tTshowste results of mixing wild-type (WT) and mutant (M) template molecules at input ratios ranging 

Sono?Toofono^ 

o^ildUe and mutant product after amplification (output ratios). These results show that TMBSP ana.ysis of melting 
curves can distinguish 1 molecule in 1 00 following 45 cycles of PCR ampl.f .cat.on. 



Precision of the assay 

[0082] Table 1 summarizes the accuracy of quantification by melting-peak analysis, ^f*^ 

discernable with reasonable accuracy and at a 100-fold difference the m.nor spec.es can st.ll be rout.ne.y detected, 

but with poorer accuracy. 

Table 1 . 



Ratios of Mutant and Wild-type alleles calculated from melting-curve analysis 




Amount of input wild-typo compared to mutant 








equal 


2X 


5X 


10X 


20X 


50X 


100X 


Mutant Copy 
number 


M/WT Ratio 
st dev 


M/WT Ratio 
st dev 


M/WT Ratio 
st dev 


M/WT Ratio 
stdev 


M/WT Ratio 
st dev 


M/WT Ratio 
st dev 


M/WT Ratio 
st dev 


to 6 


1.040 
0.041 














10 5 


1.030 
0.069 


0.544 
0.044 


0.230 
0.008 


0.108 
0.007 


0.052 
0.008 






10 4 


1.010 
0.055 


0.517 
0.009 


0.227 
0.006 


0.117 
0.010 


0.062 
0.005 


0.027 
0.011 


0.012 
0.007 


10 3 


0.943 
0.068 


. 0.503 
0.034 


0.216 
0.015 


0.104 
0.005 


0.051 
0.003 


0.034 
0.005 


0.018 . 
0.003 


102 


0967 
0.173 


0.493 
0.036 


0.207 
0.030 


0.116 
0.018 


0.058 
0.006 


0.022 
0.004 


0.011 
I 0.006 
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100831 Because of the exponential nature of PCR. small differences in reaction efficiencies will have , «*r greater 
enough to affect product quantification. 



Software 



w 



15 



20 



25 



rnnndi Current analysis software used to assess the data takes melting curve data, differentiates with respect to 

o floa For quantification wrth internal standards, the number of curves is illustratively two, and the means are known 
ooewith^ 

free The non-linear regression software can be modified to allow the user to enterthe expected melung temperatures 
of taroet and competitor and the concentration of the competitor in each sample. t „ tharnm 
moHT The le melting peak areas are used to calculate HER-2/neu target copy number. Users enter the com- 
pel copy numbe foTeach sample. The software takes the data from multiple samples and plots the log of e n 
SeZmpeTr ratio versus thefog of the competitor concentration. This plot should give a hne a s.ope of -1 w.th a 
y-intercept equal to the log of the initial target concentration. 

EXAMPLE 5 

[0087] The following experiment is conducted to determine quantification of HPV 16 using internal quantification 
standards with real-time fluorescence PCR on the LightCycleKB 



DNA/oligonucleotides 

3 o [0088] Human Papilloma vims DNA is subcloned into pBR322 . The following probes and primers are used for cloning, 
amplification, and detection. 

16HI13: 5*-GGGGATCCACTTCAGTATTGC-3' (SEQ ID NO. I); 



35 



16R19: 5 ' -GGGA ATTCC ATGGCTG ATCCTGC AGGT AC-3 (SEQ ID NO.2); 



40 16ICS: 5 ' -GATCCTGC AGGTACCGATCGGAT AGTGAGCGAGAGAT AGGTA 

GGG ATGGTTTTATGTAG-3 ' (SEQ ID NO. 3); 



45 



50 



55 



lCS P 913/640: 5'-LC640-CTACCTATCTCTCGCTCACTATCCATC-P-3' 
(SEQ ID NO.4); 

l6p913:5'-LC705-ATTACATCCCGTACCCTCTTCCCCATT-p-3' 
(SEQ ID NO.5); 

900fl6: 5 »-CCATGGCTGATCCTGC AGGT AC-3 ' (SEQ ID NO.6); 
1300rl6:5'-CCACTTCAGTATTGCCATACCC-3' (SEQ ID NO.7); 
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TCTGATATAGCATCCCCTGTtTTTTTTTCCACTA 



C AGCCTCT AC AT A AAACC-FTTC-3 ' (SEQ TD N0.8) 



Fluorescent Dyes 
Reactions 

t 0090, An artificia. system tor the detection o, initial ^P'** 

genomic DN A that had been previously cto£«to. ^^^^^XZ^ endonucLse site 
readily cloned into a pUC19 plasmid. containina plasmid DNA using a^ 

p,« mM DMA ~ ( ^/^ — - * 

quantification standard probe-binding site. Pnmers 16RI9 and ^^^J^^^ int0 a pUC i 9 plasmid. 
the directional subcloning of this amplicon into a pUC1 9 plasmid. 

Producing the Artificial IQS and HPV 16 plasmids 

[0 0 92] HPV1 6 p.asmidD N AaM0 7 co^ 

the fo-owing fina, concentrations albumin , sooCM each dNTP, and KlenTaq 

primer, 50mM Tris pH 8 3 (JSC). 4.0 CM MgM, cu|ar probes) . Thermal conditions 

2 seconds, extens.on a 78 p'^^^^t^V^ excep t for the transition between the anneaiing and extension 
each transition was set to the maximum of 20 C/secona, excspuo ethidium bromide) 

step at 10-C/second. Reactions were run on; . £™ *^ bands were excisedfrom 

gel at 80mA for 1.5 hours. The reaction products we re ^^^J^^o 42600, Beveriy, MA) according 

0.8% agarose gel, and purified as described above. - restriction endonuclease digested 

[0 093] Purified artificia. HPV 1 8 template. ^^^^^ I man acturers directions. Fo.lowing 
with Eco Rl and Bam HI (Boehnnger Mannehem B,oc ^^' S) D ™'2 Lcribed above. Purified template DNAs 
digestion. DNA was ^< *^ Biocbemicais) at 14'C over- 
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« h™i«rt fnr 5 minutes and DNA concentration was determined by spectrophotometer readings at 

(lCSp913). 

Amplifications for quantification analysis 

Ton, 0.4„ M OOOMS ,«£ Jjjnjr. ^T."^* ^ 
rescence acquisition step that was at 0.4°C/second. 



Results 



[00S7] As indicated abo.e in th. derivator, of .he del,, OX equation, ^ detecdon e»*ne, oj ^^t, "tnl 

lSl DN Ra to shows the detection efficiency of internal Quantification Standard and Artificial HPV 16 tempiate. The 

v y /» HPV 1R nNA ls s h 0 wn as detected in channel 3. As expected wun a *\nyic 

2 00] The cycle shift that occurs during the amplication whan differing initial copies of targe. D ^^^ r 
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15 



20 



25 



30 



35 



temolate concentration of 1x10*. The internal quantification standard is at initial template concentrations of 1x10= (Fig. 
1 S? 1 MS. and 1 x1 03 (Fig. 1 2c). Accordingly, the initial copies of .OS DNA range from ten fold greater ^han 
the HPV 16 (Fig 12a to ten fold less than the HPV 16 DNA (Fig. 12c). As can be seen .n F,g. 12b, where the .n.t.al 
lonv numbe oS larget and the competitor DNAs are the same, the crossing thresholds are .dent.cal. However 

Quantification standard is earlier or later, respectively, than the cycle threshold for the HPV 16 DNA. 

?oToil The cycle shift for copy number differences between the competitor DNA (ICS) and the target DNA (HPV 6) 

was p ottL as'ne change in cycle threshold, between the IQS in channel 2 and the HPV 1 6 in channel 3, versus the 

Too of the iniS. copy number of the IQS in the reaction. Fig. 13 represents the data from two separate expenmen s for 

each target, HPV 16 DNA concentration each performed in triplicate. HPV 1 6 initial template concentrates 

1 x?03 1 xl 1x105 and 1x1 0 6. Error bars are determined from the standard deviation from four mdependent react.on 

da a po nts The tandard deviations are hence a combination of intra- as we., inter-experimental variation. The . majorrty 

of he cyc.e threshold error arises from inter-experimental variation. The lines plotted represent the trendhnes to _the 

averaged delta C.T data points for each IQS and HPV 1 6 concentration. A trendline for the averaged delta C.T. data 

from a least squares analysis of the best linear fit to the points. Table 2 presents the equations for the hnear f.t to the 

roTo2i neS The analysis of delta C.T. data from internal quantification standard curves is shown in Table 2 The trendline 
Sons used to calculate target concentrations are shown with the log |T J indicated in bo.d. Amphf.cat.on efheences 
and the actual and calculated concentrations of HPV 16 target DNA are also indicated. 

Table 2. ; 



Linear best-fit trend line to the data in Figure 8. 



Trendline Equation logC Q = logE(@n) + 
logT c 


Efficiency 


Calculated Copy # 


Actual Copy # 


% Error* 


y = 0.2802X . 5.9347 
R2 = 0.989 


1.91 


0.86x10 6 


10 6 


13.96 


y = 0.2872X + 4.9235 
R2 = 0.9911 


1.94 . 


0.84x10 s 


10 s 


16.15 


y = 0.2826x + 3.9527 
R2 = 0.988 


1.92 


0.90x1 0 4 


104 


10.32 


y = 0.2944x4-2.9621 
R2 = 0.9885 


1.97 


0.92x1 0 3 


10 3 


8.36 


y = 0.289x + 1.8822 
R2 = 0.9922 


1.95 


0.76x10 2 


10 2 


23.76 
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.The abso.ute difference between ,he observed copy . and ca.cu.ated copy I ,s represented as ,he % Error .or each part.cu.ar .nR*. targe, copy *. 

T01 031 The log of the slopes of each of these lines was calculated to produce the average reaction efficiency, and 
hl°oq of the y-fmercept was used to calculate the observed target DNA concentration. The observed target concen- 
Z Z lere I e ac h o' the samp.es no greater than 24 percent varied from the estimated concentrate based on 
Umit^g dTution determination of the DNA concentration and subsequent.y the initial temp.ate copy number ,n each 

[oi04] n The use of a common set of primers to amp.ify simi.ar templates and two hybridization *> 
detec the products of those temp.ates apparent.y results in samples that have simi.ar crossmg thresholds^The apph- 
ca Son of this two color detection system to interna, quantification standards has been fac.l.tated by lhe *^n o 
an equation that uses only a simple manipu.ation of the crossing thresho.d data to produce mterna. quant.f-cat.on w.th 
a minimum dynamic range of 1 0 fold on either side of the target DNA concentration. ; „ tK , „ thor 

mTo 5 r Whife the above examples have incorporated FRET oligonucleotide probe systems, .t ,s understood that other 

pTeUem^ 

be used eliminating the need for the anchor probe. The follow.ng example uses both the FRET ol.gonucleot.de probe 
system (Sensor and Acceptor) and the single-labeled probe system (Sensor probe only). 



17 



BNSDOCID: <EP 128S314A2J_> 



EP 1 288 314 A2 



EXAMPLE 6 



[0,06, This «*i d.n,ons«.s • <» = < ^S^JiSirJSJZE! 

using aTbs,™d»nami C Modeling based Signal f ^^^^.^" S ^^Scm L 

iTe^^r,,^ 

cell types in a mixed tissue sample. 



Model Bi-allelic System 



10 



15 



20 



25 



tOIOT, Asing—depo^n^ 

^^^^^ 

;^™tS"c)^ S£ "1 FRET pair pnob. sys.ann is obs.-ed in „a F 2 enanne, 



Target Preparation 



(0,09, Tangeis^h^ype.ndn^an.alie,..^ 



ratios. 

30 Melting Curve Analysis 
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,„„., Tne naacoon n,i,«we consis.s of DMA ^ ^0^"^' SKJSSX 

Biochemicals). 

Thermodynamic Modeling Based Signal Processing (TMBSP) Algorithm 

[0111] This algorithm couples digital signal processing with a ^r™^ 3 ™^ 
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™i., t . 0 nkernal Examplesof this type of method are found in U.S. Patent Nos. 5,273,632,5,748,491 , and 5,346,306. 

^ the deconvoiution kerna. as a component o, determ.n.ng the mass 

Sec^rpriorfechnologies determine the desired quantities one at a time. Once a component of the signal 
01 11 3] ?«»nd jpnortecnnoiog the signal and determine the next component. An example of 

5 f ^7^^^^^o. ™* * - d,sc,osure determine the mass 

^^>SS!Sy^"^ Signal processing have been involved in £ aria** ; - ;- ? so, PCR 

S Vhl Fas Fourier T ansform (FFT). The basic idea is to represent the signal as a linear combination of s.nusoida 
^ SKSS^i to keep only those sinusoidal basis functions that contain reliable infom.at.on about 

20 loiim^A DFT uses a finite number of terms from a Fourier series to approximate a periodic function. The Fourier 
^rep^esenTp^c function with a reasonable amount of smoothness. Suppose ,(T) is a fluorescence me.t.ng 
signal, then the Fourier series of f(T) is 
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where the temperature is rescaled by the change of variables 



T-T mln 



Imkx 



35 T -T • 

1 max ' mm 
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The variables g(k) are the discrete Fourier coefficients of f(T). Each coefficient g(k) is calculated by computing the 
integral 



' ** 

,. H , ha per i< an efficient method of computing a set of these integrals. 

mm "o S o^ * «e number of Ihj. .™™b« computed, and some terms „e meanmgleee because 
r ' ™!S.?n the . Lai The DFT provide, a simple method o[ eliminating the ooisa torn a signal by setttng 
^Sl e^a.ad w„h the noL ecu., to zero. One Is tan I. deolde „hieh eo.tticients 

is equal to the norm of the function f(T), or 
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l/(r)|f= 



101211 Assuming the noise in the signal is small, a common and safe method of eliminating no.se .s o ,use the 
' property aTkeep enough terms of the DFT so the norms of the aotua. signal and processed s.gna. are dos >£ one 
another Specifically, one sets g(k) = 0 if g(k) I < o where o is a small tuning parameter. If K(a) .s the set of d.screte 
Fourier coefficients that have not been set to zero, then the processed s.gnal .s 



10 



15 



2 ak 



and it has the property that the processed signal is close to the actual signal since 

20 *eA'(<r) . 

is small by construction. Additionally, the periodic basis functions that only contribute ^J^^^ 
LnTare ianored Usually these basis functions oscillate rapidly and are .denuded with no.se. Th.s procedure has the 
:Zd be^that it rpplimates the actual signal with a sma.l set of data. The only data that need to be stored are 
the wave numbers in the set K(o) and the corresponding discrete Fourier coeffic.ents. 

Thermodynamic Modeling 

r0122l The fluorescence signal of a PCR product decreases or increases when the product denatures. This process 
fa ase JaSon "at can be understood using thermodynamics. The thermodynamics of phase trans.t.ons of a 

mixed materia, are based on the thermodynamics ^^f^^^^^s of subs tance i 

r0123l Consider a mixture of many substances, labeled 1,2 N. It fa* I ) is tne ^iqds nee « y 

as a function of temperature, then the Gibbs free energy of a mixture of these substances is given by 
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mixjj 

i>J 



" „h.re m , is ,he mass f r.cfions of fhe sobsfance /. Th. energy from mixing is 40^ and II accounts for changes in 

l °"l^:Z7^n7Zi a ntncion of rsmpera.ur.. .empsrafure i,se» can he rnoogn, o, as a funchon o, fioor.s- 

oeratures near a phase transition. Since temperature is a function of f. compos.ng GfJ) w.th T(f) implies that G, is a 
^rZn oiZrlLnoe. Fina.ly because G,,s a monotone function one can " & 

[0126] These observations suggest that' one can model the fluorescence of a.PCR m.xture as 



55 



where f,is the fluorescence melting signals of species /, and where f mix is the fluorescence me.ting signal of the mixture 
[01 S 2TG!ven h the fluorescence melting signa.s, f, f mi , and ignoring the fluorescence mixing terms, a good approxi- 
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nation for the mass fractions o. the substances can be found by minimizing the following objective function over ail 
choices of mj between zero and one. 
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Basis Function Approximation 

mi281 The fluorescence melting signal to be analyzed may or may nothave standard melting signal such as I, de- 
scribed above When the standard signals are missing, they must be approximated. The illustrate aPP— on 
scheme U based on the observation that the fluorescence melting signal o. the products . essentially I "ear at tenv 
^Sres greater than the me.ting transition (i.e., at the "melted phase"), and non-linear at temperatures less than the 

toon- reasonable signal, and an approximation of the remainder data is calculated. The ma.hemauca. mode, used ,s. 

rcn = p, cnii - Z, m J M J p > m J M J ^ 

and iho terms in the model are defined as the following: 

• T- temperature ■ 

fr(T) - approximated fluorescence of the melting curve 
30 . P ( T) - nonlinear function representing the fluorescence in the annealed phase 

. P4T) - linear polynomial representing the fluorescence in the melted phase _ temoerature 

. . M(T) - fraction of species j that has melted; M/T;=0 implies species, .s annealed and T ; .s the melt.ng temperature 
. m - constant mass fraction of species j in the sample vessel 

Finally all the summations in the model include N terms, with one term for each melt. 

meltmg temperature 7,a. MUM curve ia steepest in the melting transition, and ,h. widtlt at me molting 

transition wj. 
Algorithm Coupling 

r013il In the first step of the full algorithm, scaling is performed, andthe samples with no melting signal Chives-) 
are found Then parameters that specify the remainder functions are found. Finally an .terat.ve process adds one 
so remainder at a time, and minimizes the objective function 



55 



min 

in, ,<r 



,0 aime.t.n.ads.y lind the toothing parameter ,. the mass tractions ot tho kn.vm standards m, and tha molting 
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10 



15 



«:i76 of the noise in the siqnal Other methods of selecting the tolerance could be used. 

[U132 aZ chart of h full algorithm process is shown in Fig 18. The top box is the entry ,nto the a.gonthjr .and 
users soecifvThe standards and unknowns. The second box determines the scale factors of all the Signals and deter- 
mines ^S^^S^X^ are negatives. The third box signifies where the parameters of ^^^T^ 
Z deSm ned If combinations of the standards adequately mode, the unknown melt.ng c ^ Ue ^^ £ 
oLrameleTs wi! be zero In the absence of known standards t h then the approximated curves r will be used exckJS.ve£ 
SS^rJSSTLe boxes fonrn the iterative algorithm to find a., the components of the unknown^ F^t the 
Sization problem defined in this section is solved for the current set of standards and reminders. Therr .he mode 
£ cornered against the unknown, and if the fit is within a tolerance limit, the algorithm stops and reports Js resul s_U 
L fiTis no? wtt^ the tolerance limit, then the algorithm determines a new standard and repeat, the solutton of the 
minimization problem. - 
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Results 

r0134] The wild-type al.e.e (raction is estimated (as -Output") using (1) the Th . emo ^ namto f ^f^f^??Tn 
Processing (TMBSP) method, and (2) the melting peak area ratio analysis method by JJ^ 1 ^^? 
Tables 3 and 4 those outputs are compared against the actual allele fraction ,n the sample ( Input ). Ou tpu o TMBSP 
Inaysis fs obtained for all al.ele fractions regardless of the probe system used. The values agree well w. .the ^put 
values Outpul of melting peak area ratio analysis is obtained only for the FRET pair probes .n wh.ch ahele fractions 
are^groa?e than 1 oi andLor than 90°/o (Table 3). The LCDA software used for this analysis was > unab o , c .provide 
Zl Lltino oeak area ratio for the Sensor probe-only system due to the oppos.te d.rect.on of s.gnal change n the 
ZlZTZe ^V^ 16b!l7b). The LCDA software is also not able to detect allele fractions of 10% and lower, 
or' 90% and greater. 



TABLE 3. 



Data obtained using Sensor probe only 



Input fraction of wild-type allele (%) 


Output of melting peak area ratio analysis 
(%) 


Output of TMBSP analysis (%) 


2 

5 

10 

20^ 

50 

80 

90 

95 

98 




2 

5 
11 

. 20 
49 
78 
85 
92 
94 


TABLE 4. 


Data obtained using FRET pair (Sensor and Acceptor) probes 


Input fraction of wild-type allele (%) 


Output of melting peak area ratio analysis 
(%) 


Output of I MBSP analysis (%) 


2 

5 

10 

25 

33 

40 


7 
.16 
29 
33 


3 
7 

14 
24 
34 
39 
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TABLE 4. (continued) 


Data obtained using FRET pair (Sensor and Acceptor) r. 


>robes 


Input fraction of wild-type allele (%) 


Output of melting peak area ratio analysis 
(%) 


Output of TMBSP analysis (%) 

: 52 


50 
60 
66 
80 
90 
95 
98 


47 
56 
67 
78 


61 
70 
80 

92 
95 



eating a bias with the melting peak area ratio method. 
EXAMPLE 7 

[0, 361 ™s nnpk demons.™..* .h., .he gen. dosag. in a mix.ura can be qaand.isd nsing ft. Thermodynamic 
^riSJS^SSS^S^ is s.edied ,o. de,e,,dn o, Caption a.ing ,h. addi.ion », a 

■*SS25e^SSS£s=ss^ 

normalized ratio of 1 .5. 
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EXAMPLE 8 



10.. malize the dosage of the gene locus. 
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EXAMPLE 9 



[0 142] The system described in Example 7 is further simplified by the additional use of ^^^^Zn 

20 Example 7. 

EXAMPLE 10 

roi 431 This example demonstrates that the mass fraction (or molar ratio) of two or more nucleic acids in a biological 
UblTo low defection of both genes by one detection channel on the LightCycle,® 

have diHeren, melting temperatures (Tm) so that the two genes can be ^^J^^J^p^S 
2/neu probes in which the LCRed640 dye is used on the reporter probe hat has e Tm of 64 C U.8. Paten; 

Ze Tan^c i as 1 0 It is also contemplated that using similar approaches, mass fractions (or molar rat.o) of more 
SISSLaZ^It w^n the scope and spirit of the invention as described and defined in the fo.iow.ng c.a.ms. 
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SEQUENCE LISTING 

.<110> university of Utah Research^ Foundati on 
Idaho Technology 
Eyre, David 3. 
Rasmussen, Randy P. 
Cap! in, Brian E . 
Wade, Stevenson 
desilva, Deepika M. 

<120> REAL-TIME GENE QUANTIFICATION WITH INTERNAL STANDARDS 
<130> 7475-71277 

<150> US 60/316,614 
<151> 2001-08-31 

<160> 12 

<170> Patentin version 3.1 

<210> 1 
<211> 21 
<212> DNA 

<213> Artificial Sequence 

■> 

<220> 

<223> synthesized sequence for cloning. 
<400> 1 

ggggatccac ttcagtattg c 
<210> 2 

<211> 29 * 
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<212> DNA 

<213> Artificial sequence 
<220> 

<223> synthesized sequence for cloning. 

<400> 2 29 
gggaattcca tggctgatcc tgcaggtac 

<210> 3 

<211> 59 

<212> DNA 

<213> Artificial sequence 
<220> 

. <223> Synthesized sequence for cloning, 
gatcctgcag gtaccgatcg gatagtgagc gagagatagg tagggatggt tttatgtag 59 

<210> 4 

<211> 27 

<212> DNA 

<213> Artificial Sequence 



<220> 

<223> synthesized sequence for detection of internal quantification sta 
ndard. 

<220> 

<22i> misc_feature 

<222> (1) . . CI) 

<223> 5 f -LC640 Fluorescent label 



<400> 4 27 
cta.ccta.tct ctcgctcact atccatc 

<210> 5 
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<211> 27 

<212> DNA 

<213> Artificial Sequence 
<220> 

<223> synthesized sequence for detection of artificial hrv 16 sequence. 
<220> 

<221> misc_feature 

<222> CI) - • CI) 

<223> 5'-LC705 Fluorescent label 



<400> 5 ' 27 

attacatccc gtaccctctt ccccatt . 

<210> 6 

<211>. 2 2 

<212>. DNA 

<213> Artificial sequence 



<220> 

<223> synthesized sequence for cloning. 

35 <400> 6 22 

ccatggctga tcctgcaggt ac 

<210> 7 

40 <211> 22 

<212> DNA 

<213> Artificial Sequence 



<220> 

<223> Synthesized sequence for cloning. 

<400> 7 22 
ccacttcagt attgccatac cc 
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<210> 8 

<211> 59 

<212> DNA 

<213> Artificial sequence 



<220> 

<223> Synthesized sequence for detection of internal quantificati 
ndard and artificial hpv sequence. 

<220> 

<221> misc_feature . 

<222> (59).. (59) 

<223> 3'-fluoresein Fluorescent .label 



c?cgtca?ct gatatagcat cccctgtttt ' tttttccact acagcctcra cataaaacc 

<210> 9 

<211> ,23 

<212> DNA 

<213> Homo sapiens 



<220> 

<221> misc_feature 

<222> (23).. (233 

<223> 3 '-fluoresein Fluorescent label 

<400> 9 

gttcctgcat gggcggcatg aac 

<210> 10 

<211> 33 

<212> DNA 

<213> Homo sapiens 
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<220> 

<221> mi sc_feature 

<222> CD-- CI) 

<223> 5*-LC640 Fluorescent label 



<400> 10 33 
10 ggaggcccat cctcaccatc atcacactgg aag 
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<210> 


11 


<211> 


18 


<212> 


DNA 


<213> 


Homo sapiens 


<400> 11 

gcgcactggc ctcatctt 


<210> 


12 


<211> 




<212> 


DNA 


<213> 


Homo sapiens 


<400> 


12 



18 



■ - 17 
ggtcagcggc aagcaga 



Claims 

1. A method of determining mass fractions of first and second target nucleic acids present in a test sample, said 
method comprising the steps of 

(a) contacting the target nucleic acids with a fluorescent nucleic acid indicator, the indicator being configured 
to provide a signal related to the quantity of indicator hybridized to the target nucleic ac.d, the md.cator further 
configured to discriminate the target nucleic acids based on melting temperature, 

(b) illuminating the tesl sample, 

(c) monitoring fluorescent change to generate a melting curve, and f( . . rn . 

(d) using a thermodynamically based signal processing algorithm to determine the mass fract.on of the target 
so nucleic acids. 

s 

2. The method of claim 1 further comprising the steps of 

defining standard melting curves f f for each target nucleic acid having a standard melting curve, and 
55 defining a melting curve f mix of an unknown sample : 

wherein the thermodynamically based signal processing algorithm simultaneously finds a smoothing param- 
eter o. mass fraction m,of each of the nucleic acid targets present in the unknown sample, and melt.ng temperature 
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and mass fraction parameters of remainder functions, and uses a coupled optimization and iterative process to 
minimize 



10 



15 



said optimization and iterative process repeated until a sum of the mass fractions is > 1 - 6, where <= is a tolerance 
value. 

The method of claim 1 wherein the mass fraction of the target nucleic acids provides information concerning a 
deletion or a duplication in a gene. 

The method of any of claims 1-3 wherein the fluorescent nucleic acid indicator comprises a f)uorescently-labeled 
• sequence specific oligonucleotide probe. 

5. The method of claim 4 wherein the sequence specific oligonucleotide probe is selected from Ihe group consisting 
20 of a FRET pair probe system and a single-labeled oligonucleotide. 

6. The method of any of claims 1 -5 wherein the second target nucleic acid is a competitor of the first target nucleic 
acid for the fluorescent nucleic acid indicator. 

25 7 The method of any of claims 1 -6 wherein the test sample further comprises a thermostable polymerase and a pair 
TJ^S^S^ configured, for amplifying the first target nucleic acid, said method further compnsing 
the step of 

amplifying the target nucleic acid, said amplifying step occurring prior to the step of using the thermodynam- 
ic ically based signal processing algorithm. 

8. A method of quantifying a target nucleic acid present in a biological sample, said method comprising the steps of 

(a) combining in a single reaction vessel at least a portion of said sample, a thermostable polymerase^ known 
LncentratioLacompetitornu^ 

wherein said pair of oligonucleotide PCR primers is configured for amplifying a selected segment of the 
target nucleic acid and the competitor nucleic acid; wherein said competitor nucleic acid has a «"'<«"«^! on 
having a different sequence from a corresponding region of the target nucleic acd; and where.n the competitor 
nucleic acid and the target nucleic acid are amplified with essentially equal efficiency, 

said oligonucleotide probe labeled with a first fluorophore and configured to hybridize to the unique 
section of the competitor nucleic acid and the conesponding region of the target nucleic acd; 

wherein hybridization of the oligonucleotide probe to at least one of its respective complementary target 
nucleic acid and competitor nucleic acids results in a change in the magnitude of fluorescence from the fluor- 
ophore; . . „ . 

(b) amplifying the selected segment of the target and competitor nucleic acids; and 

(c) illuminating the biological sample and monitoring fluorescent change from the first fluorophore. 

9. The method of claim 8 wherein the illumination step includes melting curve analysis. 

50 10. The method of claim 9 wherein the melting curve analysis includes using a thermodynamic modeling based signal 
processing algorithm. 

11 ' further comprises a first section having a sequence identical to a corresponding 

nucleic acid and to the first region of the target nucleic acid, adjacent to the unique sect.on of the compet.tor nucle.c 
acid and adjacent to the second region of the target nucleic acid; 
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firs, fluorophore and the third fluorophore in a resonance energy transfer relauonsh.p. 
12. The method of ciaim 8 wherein step (c) inc.udes monitoring the fluorescence as a function of time, the method 
10 further comprising the steps of 

(d) creating an amplification profi.e for the target nuc.eic acid and an amp.ification profi.e for the competitor 

^Tcomraringthe amp.ification profi.e of the target nuc.eic acid with the amplification profi.e of the competitor 
ifudeteS I to determines a cyc.e shift between the competitor and target nuc.e.c ac.ds. 

: 13 The method of c.a.m 12 wherein the steps (a) through (c) are repeated using varying known concentrations of the 
' competitor nucleic acid, and the method further comprises the step of 
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